DEVELOPMENT  OF  AN  ECG  IDENTIFICATION  SYSTEM 


Masaki  Kyoso  ,  Akihiko  Uchiyama" 

'Department  of  Information  and  Computer  Sciences,  Kanagawa  Institute  of  Technology,  Kanagawa,  Japan 
2School  of  Science  and  Engineering,  Waseda  University,  Tokyo,  Japan 


Abstract  -  The  purpose  of  this  paper  is  to  present  a  new 
identification  engine  for  an  ECG  identification  system.  The 
system  identifies  subjects  based  on  a  comparison  of  a 
patient’s  ECG  with  previously  registered  ECG  feature 
parameters.  These  feature  parameters  are  sampled  from  the 
intervals  and  durations  of  the  electrocardiographic  waves 
extracted  using  characteristic  points  appearing  on  the 
waveform  of  the  second  order  derivative  and  are  identified 
using  discriminant  analysis.  Suitable  combinations  of  feature 
parameters  are  selected  and  the  accuracy  of  the  proposed 
technique  verified  using  ECGs  from  nine  normal  subjects. 
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I.  Introduction 

The  widespread  use  of  telemedicine  and  home  care 
techniques  has  increased  the  need  for  systems  to  protect 
patient  data  from  leakage,  theft  or  modification  during 
transmission.  Furthermore,  a  robust  identification  system  is 
needed  to  ensure  that  the  correct  information  is  stored  for 
each  patient.  ECG  is  a  frequently  used  method  for  monitoring 
and  checking  a  patient’s  condition.  An  identification  system 
that  utilizes  some  of  the  ECG  feature  parameters  would 
improve  the  level  of  security  and  reliability  of  the 
information  system  without  adding  any  data,  a  vital 
requirement  for  minimizing  data  volumes.  For  this  purpose, 
we  proposed  a  simpler  and  more  practical  system  [1]  than 
previous  systems  [2,3].  In  this  study,  we  improve  on  the 
identification  engine  using  multiple  discriminant  analysis  and 
demonstrate  the  performance  of  the  discrimination  process  by 
experiments  using  ECGs  from  nine  subjects. 

II.  Method 

In  order  to  perform  any  identification  process,  feature 
parameters  must  be  extracted  from  the  ECG.  In  this  study,  we 
used  the  durations  and  intervals  of  the  electrocardiographic 
waves.  Four  parameters  (P  wave  duration,  PQ  interval,  QRS 
interval  and  QT  interval)  are  appropriate  for  this  purpose, 
since  they  are  not  affected  by  the  R-R  interval  or  electrode 
conditions.  The  extraction  of  these  parameters  is  performed 
by  marking  characteristic  points  on  the  ECG  using  the  second 
order  derivative  waveform  to  emphasize  their  locations.  An 
example  of  such  marking  is  shown  in  fig.  1 . 

Discriminant  analysis  with  two  feature  parameters  is  used 
to  discriminate  between  the  registered  ECGs.  A  combination 
of  two  of  the  four  feature  parameters  are  used  as  the  variate 
for  analysis.  Since  testing  showed  that  most  covariance 
matrices  calculated  from  the  feature  parameters  of  subjects 


Fig.  1 .  Marking  at  characteristic  points  on  ECG. 

(a)  :  ECG 

(b)  :  Second  order  derivative 


are  different,  we  applied  Mahalanobis’  generalized  distance 
to  the  criterion  for  discrimination.  The  average,  variance  and 
covariance  values  are  calculated  from  the  first  sequence  of 
feature  parameter  data  to  fix  the  coefficients  in  the  equation 
of  the  Mahalanobis’  generalized  distance  for  registration 
purposes.  The  discriminant  analysis  result  is  computed  by 
selecting  the  smallest  Mahalanobis’  generalized  distance. 
Using  this  approach,  the  system  can  discriminate  between 
registered  subjects  within  a  single  beat. 

III.  Experiment 

The  experimental  system  is  divided  into  three  parts.  The 
first  is  an  ECG  measurement  and  AD  conversion  block,  the 
second  is  a  feature  parameter  extraction  block  and  the  third  is 
a  discriminant  analysis  block.  The  feature  parameter 
extraction  block  extracts  four  feature  parameters  from 
sampled  ECG  data  sequence  on  a  beat-by-beat  basis,  and 
checks  whether  the  values  are  within  the  allowable  ranges  of 
normal  electrocardiographic  waves.  Two  sequences  are  then 
selected  from  these  feature  parameter  sequences  and  applied 
to  the  final  block.  Finally,  a  beat-by-beat  discrimination 
result  is  obtained. 

The  experiment  was  performed  using  ECGs  from  nine 
normal  subjects.  Two  sequences  of  ECG  data  were  measured 


Report  Documentation  Page 


Report  Date  Report  Type  Dates  Covered  (from...  to) 

25  Oct  2001  N/A 

Title  and  Subtitle  Contract  Number 

Development  of  An  ECG  Identification  System 

Grant  Number 
Program  Element  Number 

Author(s)  Project  Number 

Task  Number 
Work  Unit  Number 

Performing  Organization  Name(s)  and  Address(es)  Performing  Organization  Report  Number 

Department  of  Information  and  Computer  Sciences  Kanagawa 
Institute  of  Technology  Kanagawa,  Japan 

Sponsoring/Monitoring  Agency  Name(s)  and  Address(es)  Sponsor/Monitor’s  Acronym(s) 

US  Army  Research,  Development  &  Standardization  Group  - 

PSC  802  Box  15  FPO  AE  09499-1500  Sponsor/Monitor’s  Report  Number(s) 

Distribution/ Availability  Statement 

Approved  for  public  release,  distribution  unlimited 

Supplementary  Notes 

Papers  from  23rd  Annual  International  Conference  of  the  IEEE  Engineering  in  Medicine  and  Biology  Society,  October 
25-28,  2001,  held  in  Istanbul,  Turkey.  See  also  ADM001351  for  entire  conference  on  cd-rom. 

Abstract 


Subject  Terms 


Number  of  Pages 

3 


for  each  subject.  The  coefficients  in  the  equation  of  a 
Mahalanobis’  generalized  distance  were  calculated  from  the 
first  ECG  data.  Discriminant  analysis  was  performed  on  the 
second  ECG  data.  Table  I  shows  the  measurement  conditions 
for  the  ECG  measurement  and  AD  conversion  block.  All 
possible  feature  parameter  combinations  were  applied  to  the 
discriminant  analysis  to  identify  the  optimum  combination. 

IV.  Results 

The  details  of  the  measured  ECGs  are  shown  in  table  II. 
The  averages  and  the  variances  of  the  feature  parameters 
from  the  measured  ECGs  are  shown  in  table  III. 

Some  examples  of  feature  parameter  distributions  are 
shown  in  fig.  2  through  fig.  4.  Fig.  2  and  fig.  3  show  results 
from  different  sequences  of  ECG  data  from  one  subject.  The 
distribution  of  the  parameters  from  a  second  subject  is  shown 
in  fig.  4. 

The  result  of  discriminant  analysis  is  shown  in  table  IV. 

Table i 

MEASUREMENT  CONDITION 


Num.  of  samples 

100,000 

Filter 

0.06Hz  HPF,  60Hz  LPF, 
50Hz  BRF 

AD  converter 

500Hz,  12bit 

Lead 

Limb  lead  11 

State  of  subject 

Rest  (sitting) 

Table  ii 

OUTLINE  OF  THE  DATA  FROM  NINE  SUBJECTS 

Subject 

Age 

Date 

ID  Num.  of  beats 

A 

21 

Oct.  19,2000 

A  1 

A  2 

218 

224 

B 

23 

Oct.  14,2000 

B  1 

B  2 

295 

291 

C 

22 

Oct.  11,2000 

C  1 

C  2 

202 

201 

D 

35 

Feb. 09, 2001 

D  1 

D  2 

170 

174 

E 

21 

Oct.  19,2000 

E  1 

E  2 

260 

275 

F 

22 

Feb. 08, 2001 

F  1 

F  2 

259 

239 

G 

22 

Oct.l  1,2000 

G  1 

G  2 

275 

255 

H 

24 

Oct.  17,2000 

H  1 

H  2 

217 

224 

I 

23 

Feb. 08, 2001 

I  1 

I  2 

227 

226 

•4ft 
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Fig.  2.  Distribution  of  characteristic  parameters,  (data  :  A_l) 
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Fig.  3.  Distribution  of  characteristic  parameters,  (data  :  A_2) 
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Fig.  4.  Distribution  of  characteristic  parameters,  (data  :  G  1 ) 


Table  iii 

AVERAGES  AND  VARIANCES  OF  THE  FEATURE  PARAMETERS 


Average 


A  1 


0.105 

0.088 

0.150 

0.369 

1.94E-5 

8.06E-7 

1.78E-5 

1.90E-5 


B  1 


0.094 

0.097 

0.142 

0.335 

6.55E-5 

1.10E-6 

4.47E-5 

2.45E-5 


C  1 


0.116 

0.095 

0.167 

0.329 

1.66E-4 

1.44E-6 

4.18E-5 

1.87E-5 


D  1 


0.088 

0.096 

0.137 

0.426 

2.48E-5 

1.01E-6 

1.87E-5 

6.27E-6 


E  1 


0.104 

0.100 

0.157 

0.354 

5.07E-5 

9.42E-7 

3.32E-5 

1.62E-5 


F  1 


0.101 

0.093 

0.147 

0.375 

7.58E-5 

1.24E-6 

7.08E-5 

2.42E-5 


G  1 


0.099 

0.090 

0.178 

0.331 

4.26E-5 

2.18E-6 

1.45E-5 

5.16E-5 


H  1 


0.105 

0.109 

0.142 

0.378 

7.15E-6 

9.89E-7 

7.44E-6 

1.27E-5 


I  1 


0.097 

0.105 

0.153 

0.337 

1.93E-4 

2.35E-6 

8.45E-5 

9.08E-5 


Variance 


The  accuracy  value  quoted  in  the  table  is  calculated  from  the 
number  of  beats  correctly  discriminated  and  the  total  number 
of  beats. 

V.  Discussions 

The  statistical  parameters  listed  in  Table  III  were  used  as 
coefficients  in  the  calculation  of  the  Mahalanobis’ 
generalized  distance  it  is  demonstrated  that  these  parameters 
are  suitable  for  discrimination  purposes  in  fig.  2  through  fig. 
4.  Fig.  2  and  fig.  3,  which  show  the  distribution  of  feature 
parameters  from  the  same  subject,  show  a  high  degree  of 
correlation.  On  the  contrary,  fig.  3  and  fig. 4,  which  show  the 
distribution  of  feature  parameters  from  different  subjects,  are 
significantly  different. 

Table  IV  shows  the  system  performance  for  every  possible 
pair  of  feature  parameters.  It  can  be  seen  that  there  is  a  large 
variation  in  achievable  accuracy  with  the  different 
combinations.  In  particular,  use  of  the  P  wave  duration 
results  in  a  low  discriminatory  accuracy.  This  may  be  due  to 
the  small  differences  between  subjects  and  from  ambiguities 
in  the  placement  of  the  characteristic  points.  The  highest 
accuracy  was  achieved  using  the  combination  of  QRS 
interval  and  QT  interval  parameters.  All  results  with  the 
exception  of  the  C_2  data  exceed  ninety  percent  accuracy. 
Since  the  combination  of  QRS  interval  and  PQ  interval 
parameters  also  provided  accurate  results,  it  is  expected  that 
the  use  of  a  combination  of  these  three  feature  parameters 
will  improve  the  accuracy  of  the  technique. 

The  ratio  of  the  two  amplitudes  of  the  electrocardiographic 
waves  is  another  possible  parameter  for  use  in  this  system, 
since  it  is  not  affected  by  variations  in  the  amplitude  of  the 
ECG  due  to  uncertain  electrode  conditions.  The  extension  of 
the  system  to  include  this  feature  parameter  may  suppress 
any  deterioration  in  performance  that  may  occur  when  the 
number  of  subjects  in  the  database  is  large. 


VI.  Conclusion 

In  this  paper,  we  presented  a  new  discrimination  process 
for  an  ECG  identification  system  and  demonstrated  the 
accuracy  of  the  proposed  technique  using  the  ECGs  from 
nine  normal  subjects.  The  results  suggest  that  the  system  can 
discriminate  between  registered  ECGs  with  an  acceptable 
level  of  accuracy.  It  was  also  suggested  that  the  addition  of 
extra  feature  parameters  to  the  discrimination  system  would 
further  enhance  the  performance  of  the  system. 
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Table  iv 

RESULTS  OF  DISCRIMINANT  ANALYSIS 


Data 

Accuracy  [%] 

ID 

P-QRS 

P-PQ 

P-QT 

QRS-PQ 

QRS-QT 

PQ-QT 

A_2 

80.8 

31.3 

48.7 

99.1 

99.1 

48.2 

B_2 

73.2 

38.5 

4.1 

89.0 

90.4 

46.4 

C_2 

84.0 

90.1 

81.0 

95.0 

76.1 

78.1 

D_2 

60.3 

56.3 

99.4 

44.8 

99.4 

99.4 

E_2 

82.9 

17.5 

45.8 

93.8 

90.2 

47.6 

F_2 

63.6 

7.5 

86.6 

67.4 

99.2 

87.0 

G_2 

64.3 

94.9 

38.4 

98.4 

96.5 

96.9 

H_2 

80.4 

84.4 

29.5 

79.0 

97.3 

24.1 

I_2 

99.6 

88.5 

83.2 

99.6 

99.6 

84.5 

